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Stiffness and hysteretic energy losses are calculated for a magnetically levitating system composed 
of a type-II superconductor and a permanent magnet when a small vibration is produced in the 
system. We consider a cylindrically symmetric configuration with only vertical movements and 
calculate the current profiles under the assumption of the critical state model. The calculations, 
based on magnetic energy minimization, take into account the demagnetization fields inside the 
superconductor and the actual shape of the applied field. The dependence of stiffness and hysteretic 
energy losses upon the different important parameters of the system such as the superconductor 
aspect ratio, the relative size of the superconductor-permanent magnet, and the critical current of 
the superconductor are all systematically studied. Finally, in view of the results, we provide some 
trends on how a system such as the one studied here could be designed in order to optimize both 
the stiffness and the hysteretic losses. 



I. INTRODUCTION 

Superconducting levitation results from the interac- 
tion between an applied field and the currents in a 
superconductor [Q. Macroscopically, a type-II supercon- 
ductor with strong pinning can be studied using the criti- 
cal state modelQ]. Within this framework, currents with 
constant density are induced in the superconductor due 
to a variation of the applied field it feels, either for a 
change in the value of a uniform field or for a displace- 
ment of the superconductor in a non-uniform field. A rel- 
ative movement between a superconductor and a perma- 
nent magnet will thus result in a magnetic force between 
the induced currents in the superconductor and the field 
of the magnet. If this force compensates the weight of 
the levitated object and the conditions for stability hold, 
the object will stably and passively levitate||. 

Earnshaw proved that only diamagnetic materials can 
passively and stably levitate [Q. Although type-I super- 
conductors can levitate (actually, they can produce large 
levitation forces because of their strong diamagnetism) , 
the use of type-II superconductors has become more use- 
ful because of their rigidity. The rigidity arises from the 
pinning of the flux lines and allows a continuous range 
of passive equilibrium positions^]. Yet, the use of the 
type-II superconductors results in an energetic loss due 
mainly to the magnetic hysteresis^. 

A superconducting bearing is a system composed, in 
its simplest form, of a superconductor and a permanent 
magnet. The rigidity and the losses are important factors 
in the design of a bearing. Actually, complicated types of 
bearings have been designed in order to improve rigidity 
and losses of the whole system. An excellent review on 
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superconducting bearings can be found in 

Measurements of stiffness j|, ^ and losses [ 
itating superconductor-permanent magnet bearings have 
been recently published. From the theoretical point 
of view, the study of stability and losses in high T c - 
permanent magnet bearings is a complicated task, mainly 
due to the trouble of finding the current penetration in- 
side the superconductor when the applied field is not ho- 
mogeneous. Several models have been presented. In |l^ ] 
and jl3| we proposed approximate analytical expressions 
for stiffness and losses, assuming both a small supercon- 
ductor and approximate demagnetizing factors. In Ref. 
H the stiffness of a bearing with a thin film supercon- 
ductor was calculated, considering the exact form of the 
applied field and assuming a fully penetrated thin super- 
conductor. Based on the interaction between dipole im- 
ages, Kordyuk Jl4j calculated the force and energy losses 
on a permanent magnet over a flat ideally hard super- 
conductor. Using the same model, Hull and Cansiz(|] 
calculated the stiffness of such a system. 

In Refs. |l5|, [l6| a model to calculate the current pen- 
etration of currents inside a finite cylinder, taking into 
account both the demagnetization factors and the radial 
components (and inhomogeneity) of the applied field, was 
presented. The model, based on the critical state and on 
the minimization of magnetic energy, was used for cal- 
culating the current penetration and levitation forces in 
the case of a cylindrically symmetric bearing system [[l6| . 
In this work, we will extend the application of the model 
to calculate, not only the current penetration and the 
levitation force, but the rigidity and the hysteretic en- 
ergy losses resulting from the reaction of the system to 
an eventual vibration. We will study the dependence of 
these properties upon both intrinsic parameters of the 
superconductor (its critical current) and geometrical pa- 
rameters (aspect ratio of the superconductor and relative 
size between superconductor and permanent magnet). 

This paper is structured as follows. In Section II we re- 
view the calculation model used, focusing on its applica- 



2 



tion to the present case of vertical vibrations. In Section 
III we define the stiffness and the hysteretic losses and 
expound some general trends of these parameters. Re- 
sults for the case of cylindrical symmetry are presented in 
Section IV. Section V is devoted to the discussion of the 
results and their use to optimize a real system. Finally, 
the conclusions can be found in Section VI. 



II. MODELIZATION 

We consider a zero-field cooled type-II cylindrical su- 
perconductor (SC) of radius R and length L located, ini- 
tially, very far from a coaxial cylindrical permanent mag- 
net with uniform magnetization in the direction of its axis 
(PM). The PM has radius a, length 6, and magnetization 
Mpm- A scheme of the system is shown in Fig. |l|. We 
shall use cylindrical coordinates with the origin located 
at the top face of the PM. We consider that the PM is 
fixed and that the SC, initially above it, can be moved 
towards or away from it. We define the distance between 
the PM and the SC, d, as the ^-height of the bottom face 
of the SC. The general process of descending the SC from 
a very distant position to a minimum distance, d m i n , and 
then ascending it to far from the PM, d — > oo, is called 
the major loop. 

We want to simulate the response of the SC after a 
small vibration of the system, produced at some point 
of the major loop. A vibration is regarded as a change 
in the actual direction of the movement and the subse- 
quent return to the original position. This movement is 
referred to as a minor loop (see Fig. Q). In this work, we 
only consider cylindrically symmetric situations. That 
is, we only consider movements of the SC along the axial 
direction, during both major and minor loops. 

The response of a finite cylindrical SC in the critical 
state and after a movement in the presence of a non- 
uniform (but cylindrically symmetric) applied field has 
been studied in ]l6[ | using a numerical procedure de- 
scribed in The model includes the following as- 
sumptions, besides the cylindrical symmetry: a) no equi- 
librium magnetization for the SC is considered; b) the 
PM is not affected by the presence of the SC; and c) 
the superconductor is assumed to be in the critical state. 
(The force creep due to thermal demagnetization has 
been recently studied by Qin et al |l7| ) . We calculated in 
p6[ the levitation force during a major loop and studied 
the principal characteristics of that force and its depen- 
dence on important parameters of the system. We now 
briefly review the main steps in the calculation of the 
desired magnitudes. After a variation of the position of 
the superconductor, we calculate the currents induced 
in the superconductor by imposing that these currents 
tend to minimize the magnetic energy of the supercon- 
ductor. When a variation of the direction of movement 
is produced, currents in reverse direction begin to en- 
ter the superconductor. In this latter case, we calculate 
the response of the SC using the typical procedure of the 



critical state 18 keeping the already set currents and cal- 
culating the penetration of new reverse ones in the same 
way as for the initial stage but considering the variation 
of the field. 

For a given current distribution inside the supercon- 
ductor, the magnetic force between the currents and the 
external applied field in the case of cylindrical symmetry 
has only a z-component given by 



F, 
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where V = wR 2 L is the volume of the SC and Hp is the 
radial component of the applied field created by the PM. 

Although in the described procedure we could imple- 
ment a dependence of the critical current with the inter- 
nal field, in this work we use a constant critical current 
J c , for the sake of simplicity. In that case, Eq. [l] becomes 



(2) 



F z =2np J c / H?(p,z)pdp 



being Q the current penetrated region. Q depends not 
only on the value of J c but on the geometry of the sam- 
ple, because of the demagnetizing fields. If the supercon- 
ductor becomes fully penetrated by currents, Q, = V and 
the vertical force is proportional to J c . 

Unless the contrary is explicitly written, we will con- 
sider the following parameters: a = b = R = L = 0.01m, 
J c = 2.81 10 7 A/m 2 , and M PM = 7.95 10 5 A/m. The geo- 
metrical values have been chosen among the typical ones 
in levitation experiments. The value for Mpm is such 
that pqMpm =1 T. J c is chosen to equal H /R, being H 
the magnetic field at the origin of coordinates (the center 
of the top face of the PM). Using these parameters, the 
field Hq can be considered as a typical field of penetration 
of the superconductor since when the superconductor is 
at zero distance (so the applied field in the center of its 
lower face is Hq) the superconductor is considerably pen- 
etrated by currents (although not completely penetrated 
because the sample is finite and H p = J C R is the pene- 
tration field for an infinite sample in a uniform applied 
field). We will see below why this is an appropriate pa- 
rameter. The use of non-normalized variables allows the 
easier comparison with experimental values. However, 
an adequate normalization could reduce the number of 
parameters of the system (see Ref. |l2|). In this work, 
we will use non-normalized variables for the calculations 
but will discuss the results in relation to important di- 
mensionless parameters. 

In Fig. U we show calculated results illustrating the 
general shape of the major loop of the force-distance 
curve together with some minor loops at several start- 
ing distances. In the inset, a detail of one minor loop 
is also showed. Since the penetration of currents in the 
critical state is a hysteretic phenomenon, so the force is, 
in both the major and the minor loops. In Fig. || we 
show calculated current penetration profiles correspond- 
ing to several points of one minor loop (we have chosen 
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the loop in the inset of Fig. [I]; the plotted points are 
marked there). We observe the hysteretic process of cur- 
rent penetration causing the hysteresis in the force during 
a minor loop. As we shall see in detail in the next sections 
one can analyze the stability and the hysteretic losses of 
the system from the shape of these minor loops. 



III. DEFINITIONS AND GENERAL 
CONSIDERATIONS 

A. Restoring force and stiffness 

A useful mechanical system, including those with lev- 
itating components, is often required to be stable in all 
directions. Any variation of the position of the working 
point of a rigid solid must result in a restoring force that 
tend to return the system to its original place. Consider- 
ing small displacements from the working position, one 
could define the stiffness as the first order spring constant 
of the restoring force. Thus, we define a stiffness matrix 
whose elements are 



Kaf3 



(3) 



The physical meaning of the a/3-element is the following: 
it is the first order spring constant of the restoring a- 
component of the force, after a variation of the working 
point ro in the /3-direction. The minus sing is set by 
convenience, as discussed below. 

In this work we consider a cylindrical system in which 
the working position and the vibrations of the system will 
always maintain the cylindrical symmetry. The force af- 
ter one of such variations will have only axial component, 
which results in K pz — ng z = 0. In that case, thus, only 
the component k zz should be accounted for. From Eq. 
H, we have that 



K zz (zi) 



OF. 



I ll I 



dz 



(4) 



The upper index ml indicates that the force that should 
be considered is the one within the minor loop. This is 
because, in the case of type-II superconductors, the force 
is hysteretic, so the force after the start of a vibration is 
the one inside a minor loop. 

Although in this work we only calculate the k zz compo- 
nent, it is interesting to discuss some general properties of 
the stiffness matrix in the case of a cylindrically symmet- 
ric working point. Since there is no angular dependence, 
any displacement of the superconductor in the angular di- 
rection should not modify the current distribution so the 
restoring force is zero. Thus, n p e = nee = &z6 = 0. More- 
over, by symmetry, any displacement in the radial direc- 
tion will tend to give a restoring force with no-angular 
component, resulting in k,q p = 0. In that latter case there 
will appear both a restoring force in the radial direction 
(and stiffness n pp ) which is the lateral restoring force, 



and also a vertical force (with stiffness n zp ) which will 
tend to modify the height of the working point. 

The stiffness coefficients can have different values for 
the same working point depending on how that working 
point has been achieved. This is due to the hysteresis in 
the behavior of the levitation force in type-II supercon- 
ductors. Although the stiffness is only calculated from 
the force during a minor loop, that force depends on the 
previous magnetic history of the superconductor. In the 
case considered here of vertical movements, we expect to 
obtain different value for the k zz (z%) for a given Z\ when 
considering a minor loop produced at a height Z\ during 
the general descending stage or when considering a minor 
loop produced in the same z\ but in the major ascending 
stage. 

Finally, the minus sign in Eq. || is set for convenience, 
since in this case the condition for vertical stability is 
n zz > 0. If, at some equilibrium position, an increment in 
the vertical distance is produced (dz > 0), the levitation 
force should decrease (dF z < 0; the weight of the object 
would tend to return it to the equilibrium distance). On 
the other hand, if the levitated object is pushed down 
(dz < 0), the levitation force should increase (dF z > 0). 
In both cases k zz > 0. 



B. Hysteresis and energy losses 

The hysteretic energy losses can be defined as the work 
done by the system during a complete minor loop. They 
can alternatively be defined as the energy needed to force 
a levitating object to follow a minor loop and return to 
the original position. The only difference is the sign of 
the two expressions. To avoid this ambiguity, we define 
the hysteretic energy losses as 



E 



Fdr 



ml 



(5) 



where the subindex ml indicates that the integration 
should be done over the minor loop. 

This general expression can be simplified when consid- 
ering that the cylindrical symmetry is preserved through 
the minor loop. In that case, F p = Fe = 0. If we con- 
sider that the minor loop runs between Z\ and z-i (with 
z\ < Z2) and name E zz the resulting hysteretic energy 
losses, Eq. |^ is simplified into 



E zz ( Zl ) 



\F}-F}\dz, 



(6) 



where Fj and Fj refer, respectively, to the force during 
the ascending and descending movements of the super- 
conductor during the considered minor loop. 

The point Z\ at which E zz is calculated is the starting 
point of the vibration only if produced during the major 
descending stage. If the minor loop is produced during 
the major ascending stage, then z\ is the reversing point 
of the vibration. The reason is that, differently from 
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the stiffness coefficient k zz , the energy losses at a given 
point do not depend (in the critical state framework) on 
the previous magnetic history of the superconductor. In 
the case considered here of vertical movements, E zz (z\) 
does not depend on whether z± is the one in the gen- 
eral descending stage or in the ascending one. This fact 
can be understood from Eq. ||. Since in the calculation 
of E zz one must consider the difference between forces 
during the minor loop, the eventual contribution of the 
frozen currents to the force is eliminated from the equa- 
tion, resulting in the non-hysteretic behavior of E zz . 



IV. RESULTS FOR n zz AND E zz . 

In this section we present calculations of the vertical 
stability and the hysteretic energy losses of a cylindrically 
levitated system. We systematically study the depen- 
dence of stiffness k zz and losses E zz on different relevant 
parameters of the system such as the critical current of 
the SC, the aspect ratio of the SC, and the relative di- 
mensions of the SC and the PM. We consider in all cases 
minor loops started at a fixed distance (z\ = 0.005m) 
and amplitude of the vibration (Az = 0.0025m). 

Both coefficients k zz and E zz have some general fea- 
tures independently of the parameters of the system (ex- 
cept in some extreme limits). During the initial de- 
scending, both k zz and E zz increase as z\ decreases 
owing to the increasing of the field gradient and the 
larger penetration As stated above, n zz is hysteretic 
whereas E zz is not. Another general feature is that, in 
the returning stage, the stiffness becomes negative for a 
large enough z\S. This fact, already pointed out in Jl2| , 
is produced because for large enough z\S the variation 
of the applied field is small and the currents hardly enter 
inside the sample. The minor loops have then to follow 
the major one, resulting in a negative stiffness (positive 
initial slope). Actually, depending of the shape of the 
PM, it may also exists a non stable region in the major 
descending loopJl2[. 



A. J c dependence 

In Fig. | we show the calculated k zz (Fig. ^-a) and E zz 
(Fig. ^-b) as a function of the critical current J c , for three 
values of the aspect ratio L/R of the SC, corresponding 
respectively to a thin film sample (L/R= 1/5), a regular 
size sample (L/R — 1), and a large sample (L/R — 5). In 
these calculations we have fixed the PM and the value of 
R (maintaining so the value a/R). We observe that the 
stiffness increases with increasing J c , whereas E zz has a 
maximum and tends to zero in both limits for low- J c and 
high-J c . 

Results for the stiffness can be understood from a sim- 
plified model that considers a small PM and a very thin 
superconductor [p| p0| . Within this model, n zz has two 



contributions: 



Kzz(zi) = -fM>TT R 2 L 



-Xo 



dH z 
dz 



+ M z (z x ) 



d 2 H z 
dz 2 



The first term in the Eq. [7] is independent of J c and is 
related to the initial susceptibility of the SC, xo = 

. The second one depends upon the magnetization of 
the SC and, thus, upon J c . Riise et al |5| reported that 
for a typical experimental setup the second term is much 
less important that the first one, so, for their very thin 
film, the stiffness was independent of J c . 

In the present general case, there are also two terms de- 
termining the stiffness. Consider a minor loop produced 
during the descending of the superconductor. The force, 
just at the starting point of the minor loop, z±, is 



Fi r ( Zl ) = -27T/i I Jc I 



Jb(z',zi) 



iJ>,z>dpdz', 



(8) 

where b(z',zi) describes the current penetration profile. 
If at this point a minor loop starts, these already induced 
current are kept frozen (the superscript fr in Eq. || in- 
dicates this) and new reversal currents begin to enter. 
Then, the force can be calculated as 



F z (d) = Ff r (zi) + F:™(d), 



(9) 



where Fl lew is the force due to the new induced currents: 

rL r R 

i(z' ,d) 

(10) 

being a(z' , d) the current profile of the new induced cur- 
rents. 

The coefficient k zz is calculated from Eq. || as 



FT w (d) = +27T/1Q I J c I 2 f L [ H a p (p, z')pdpdz', 

JO Ja(z'.d) 



Kzz(zi) 



dF/ r dFT 



dd dd 



(11) 



d >z\ 



Although the currents induced prior to the minor loop 
are kept frozen, the first term in Eq. [ll] is different from 
zero because the applied field involved in the integral 
changes when varying z. This term corresponds to the 
M z (zi) d g ^/ term in Eq. ]?]. It depends on J c . On the 
other hand, the second term in Eq. ^ depends on the 
initial slope of the magnetization of the superconducting 
sample. It is the equivalent to the non-hysteretic term 
in Eq. [7] Within the critical state model, this term does 
not depend on J c . 

Considering a minor loop started at a given working 
point, increasing J c increases the magnetization of the 
SC sample. The first term in Eq. [ll] increases whereas 
the second one does not. The result is an increase of the 
stiffness, as can be seen in Fig. ||-a. 

In the large J c limit the force is almost non-hysteretic. 
Therefore, in this limit, the minor loop runs over the 
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major one and the stiffness can be approximately calcu- 
lated directly from the slope of the major loop at a given 
distance. 

The results for E zz can be understood by looking at 
Fig. m where we have plotted several minor loops corre- 
sponding to the cases L/R = 0.2 (left), L/R = 1 (mid- 
dle), and L/R = 5 (right), for different J c values. (All 
the discussion below holds for the three figures, although 
the complete discussion is best seen in Fig. |5^a) . For the 
case that J c is small (actually, when it is small enough 
to produce a full penetration of the sample for a small 
variation in d) as soon as the minor loop is started, the 
sample is fully penetrated and as soon as the minor loop 
is reversed, the sample is again fully penetrated by cur- 
rents circulating in the opposite direction. The result is 
an almost symmetric (with respect to the zero-force) be- 
havior of the force during the minor loop. In that case, 
the levitation force is roughly proportional to the critical 
current (Eq. ||) and, from Eq. |[ we can see that the 
energy losses are 



E zz ~ Anfi | J c 



g(h) | dh, 



where g{h) is defined as 



9(h) = 



h+L 



H a Ap, z)pdpdz. 



(12) 



(13) 



For low J c , thus, E zz tends to zero, proportionally to J c . 
As J c increases, the full penetration is obtained only af- 
ter a significant portion of the minor loop has been gone 
through. The minor loop becomes less and less symmet- 
ric, but at the same time, the value of the force increases. 
The resulting area of the minor loop increases for a range 
of not very large J c 's. When J c goes on increasing, no 
full penetration of reverse currents is achieved. Actually, 
for large J c 's, the superconductor is hardly penetrated, 
either for direct currents or for reverse ones, so the losses 
are being reduced even though the value of the force in- 
creases. The resulting E zz is thus tending to zero for 
large J c , 

In Fig. ^-b we observe that the maximum in the E zz 
coefficient is displaced to lower J c when the SC is larger. 
This is understood again in terms of the penetration. For 
longer superconductors (Fig. |^-b,c), the demagnetization 
fields are lower and, as a consequence, the running dis- 
tance necessary to fully penetrate the SC is larger. The 
described crossover behavior of producing less hystere- 
sis because of less penetration is attained, thus, at lower 
J c 's. A more detailed descri ption of the dependence on 
L /R will be given in Section [V A . 



B. L/R dependence 

In Fig. 6 we plot the dependence of k zz (Fig. ^|-a) and 
E zz (Fig. ^|-b) on L/R for different values of the critical 
current. In order to study the dependence only on L/R 



we have fixed the PM, the value of the R (maintaining 
so the relation a/R), and also the starting point (z\ — 
0.005m) and the amplitude (Az = 0.0025m) of the minor 
loops. 

We observe in Fig. [j|-a that the stiffness significatively 
depends on L/R for low-J c but hardly for high-J c . The 
coefficient E zz has a singular behavior. As seen in Fig. ^ 
b, for low critical current, E zz increases for low L/R, and 
saturates for long L/R. On the other hand, for high J c , 
it decreases all the time (within the L/R range shown), 
even for small L/R, saturating for long samples. 

The argument exposed in the previous section explains 
also the behavior of k zz in the following way. The term in 
the stiffness depending on the magnetization is more im- 
portant for high-J c and, in this case, the magnetization 
slightly depends on L/R since almost no current pene- 
trates the SC specially far from the PM. So, increasing 
L does not produces significative changes (except when 
the SC is very thin). On the other hand, when J c is low, 
the sample is significatively penetrated by currents and 
the relation L/R plays an important role. 

The behavior of E zz can be more clearly understood 
by looking at Fig. ^-a and Fig. ^-b, where we have plot- 
ted the minor loops for different lengths and for a case of 
low critical current (Fig. ^-a) and a case of higher criti- 
cal current (Fig. 0-b) . When the current is low the force 
during a minor loop is al most symmetric with respect to 
zero-force value (Section IV A). Even in this low-current 
limit, when the superconductor is long, there are regions 
free of current, so the behavior departs from the sym- 
metrical one. The departure, however, is not very strong 
since the applied field in the regions where no fully pen- 
etration is achieved is weak. In Fig. ^-a we observe this 
behavior and see that E zz tends to increase as L/R de- 
creases. The appearance of a maximum in the E zz versus 
L/R curve for this case has the same origin as discussed 
in the previous section (a similar behavior of that shown 
in Fig. ^| could be observed in this case, but as function 
of L/R; the difference is that now, when L/R is large, 
the minor loops tend to saturate to a constant shape in- 
stead of changing indefinitely) . When the critical current 
is high, only when L/R is small the penetration is com- 
plete because the demagnetization effects are important 
in that case. When increasing L/R, the minor loops get 
less and less hysteretic until their shape (and thus, their 
area) is saturated (see Fig. 0-b). Actually, even in the 
case of high- J c for thin enough samples the area of the 
minor loop should tend to zero, producing also a maxi- 
mum in the force. In the scale used in Fig. ^-b, this is 
not seen. 

Since both k zz and E zz saturate for long superconduc- 
tors, their value per unit volume decrease as 1/L for long 
enough superconductors. 
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C. a/R dependence 

We study in this section the dependence of n zz and 
E zz upon the ratio between the radii of the PM and the 
SC, a/R. We have maintained the value of R in order 
to not change the relation L/R. The value of b is also 
maintained and so is Mpm- We note, however, that when 
changing a, the applied field will change. In Fig. ||-a and 
Fig. |^-b we show, respectively, the calculated k zz and 
E zz for three cases of J c corresponding to 0.1, 1, and 10 
times the relation H /R. 

In Fig. ^|-a we observe that there is a region of max- 
imum stability (maximum value of stiffness) when the 
superconductor and the permanent magnet have similar 
radius (in Ref. |l6| it was shown that also the magnitude 
of the force is larger when a ~ R). E zz shows a maxi- 
mum at a ~ R and tends to zero for both a « R and 
a >> R limits. 

The results can be explained as follows. Given a mag- 
netization Mpm of the permanent magnet, when the 
magnet has a radius a — ► the value of the applied 
field tends to zero and also does the levitation force. Both 
stiffness and energy losses should be zero in that limit. 
On the other hand, when the PM is much larger than 
the SC, in the region occupied by the superconductor 
the field is approximately homogeneous, resulting also in 
an almost null levitation force and, as a consequence, a 
very low stiffness and losses. In the intermediate regions, 
there should be, at least, one maximum. The fact that 
the maximum corresponds to a ~ R is related to the fact 
that, at a given distance from the PM, the applied field is 
most inhomogeneous at radial distances in the range of a. 
Therefore, SC's of radius R ~ a will feel larger variations 
in the applied field, resulting in a larger penetration of 
currents, which means a larger rigidity and, at the same 
time, larger losses, because of the larger hysteresis. 

V. DISCUSSION 

The levitation force of type-II superconductors is the 
result of the interaction between a distribution of induced 
critical currents and the external applied field. Stiffness 
and losses can be interpreted in terms of the current pen- 
etration profiles. In a qualitative way, the key parameter 
that accounts for the behavior of a bearing system is 
the relation between the field of penetration of the SC 
and the applied field of the PM. Other parameters are 
related with this. The penetration field for the supercon- 
ductor depends on the critical current as well as on its 
dimensions and aspect ratio. (In this work we have only 
worked with cylinders, but in general, the penetration 
field also depends on the geometry of the superconduct- 
ing sample). Actually, the larger the critical current is, 
the larger the penetration field is (proportionally) and 
the smaller L/R is, the larger that penetration field is 
(not in a proportional way H)). The applied field felt 
by the superconductor depends on the magnetization of 



the PM (proportionally, if considered uniform and con- 
stant), on the dimensions of the PM and on the distance 
between them. 

The relation between these two fields determines the 
current penetration profiles in the following way: if the 
applied field felt by the SC is much smaller than its pene- 
tration field, current will hardly penetrate the SC, which 
will yield in almost non-hysteretic behavior and large 
forces. This condition can be achieved by several rea- 
sons (or by a combination of them), namely, high J c , 
large L/R, low Mpm, small a (in relation to R), and/or 
small b. On the other hand, if the applied field felt by the 
SC is much larger that its penetration field, currents will 
penetrate a considerable amount of SC material after a 
small field variation (after a small movement of the SC). 
This will result in a highly hysteretic behavior. This con- 
dition will be fulfilled for low J c , small L/R, high Mpm, 
and/or a ~ R. 

For the quantitative description of the stiffness and 
losses one should also take into account the actual value 
of the force. For example, even when there is large hys- 
teresis, the area of a minor loop can be small if the 
applied field (and the force) is small (Fig. ||-a). Also, 
large hysteresis (almost symmetrical minor loops) does 
not necessarily means a large stiffness. We have seen that 
large J c 's increases the initial slope of a minor loop (in 
absolute value) although the hysteresis is reduced (Fig. 
|and Fig. |). 

VI. CONCLUSIONS 

Using a model based on the critical-state approxima- 
tion, we have been able to calculate the response of a 
levitating type-II superconductor in the presence on the 
non-uniform applied field created by a permanent mag- 
net. We have focused our analysis on vertical vibrations 
that can occur at a given working point. 

We have seen that when the radius of the PM and SC 
are similar, the stability of the system is maximized, al- 
though the losses are also the largest. As to the length of 
the superconductor, a very long sample can have a signif- 
icant amount of material that contribute neither to the 
force, nor to the stiffness, nor to the hysteretic energy 
losses, resulting therefore in a useless waste of SC mate- 
rial. By increasing the critical current of the supercon- 
ductor, the stability of the system is increased, whereas 
the losses could be reduced, if this critical current is high 
enough. 

Optimization of a cylindrical bearing system should 
take into account all these results. Depending on the 
needs and the restrictions for a particular design, one 
should choose a permanent magnet and a superconductor 
that fit the needs. In particular, very high critical current 
could produce large stability with low losses. However if 
a thin film SC is used, even with such high critical current 
the sample could be easily fully penetrated and the losses 
would be important. If large levitation forces are needed, 
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a long superconductor with high critical current would be 
the optimum solution for the reduction of losses, without 
considerably reducing the stability. 

In its current form, the present model is unable to de- 
scribe some important phenomena of levitation, such lat- 
eral displacements (and lateral stiffness and hysteretic 
losses), movements at high frequencies (such that the 
critical state does not have time to be established in- 
side the SC), force creep, or even finite geometries other 
than cylindrical. However, the same methodology could 
be eventually applied to these cases. If by some means, 
a procedure to obtain the current profiles inside the su- 
perconductor is found, the procedure described would be 



valid for a phenomenological description of the levitation 
force. In particular, the relation between a typical pene- 
tration field of the superconductor and a typical applied 
field is expected to be the key parameter to explain the 
experimental data, even in situations different from the 
cylindrical symmetry described here. 
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FIG. 1: Geometry and dimensions of the considered 
superconductor-permanent magnet system (left). Sketch of 
a minor loop running between z\ and £2 (right). 



FIG. 2: Vertical levitation force for a complete major loop and 
for several minor loops of amplitude Az=0.005m produced 
every 0.005m. Other parameters of the calculated case are 
in the text. In the inset there is a zoom of one of the minor 
loops. The currentpenetration profiles in the marked points 
are shown in Fig. pi 



FIG. 3: Current penetration profiles for a semi plane of con- 
stant angle (the 2-axis is at left in every figure) at several 
distances during a minor loop. The distances are the ones 
marked in Fig. ^. Currents induced during the descend of 
the minor loop (increase of the applied field) are drawn in 
dark gray. In light gray there are the reverse currents induced 
during the ascent of the superconductor. 
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FIG. 4: (a) Vertical stiffness k zz and (b) hysteretic energy 
losses E zz as a function of the critical current of the su- 
perconductor. The PM is fixed and the aspect ratio of the 
superconductor is changed by changing the value of L such 
that L/R=0.2 (dotted line), L/R=l (solid line), and L/R=5 
(dashed line). Other parameters are written in the text. 



FIG. 5: Force versus distance minor loops started at 
zi=0.005m and of amplitude Az=0. 0025m for the cases (a) 
L/R = 0.2, (b) L/R = 1, and (c) L/R = 5. The values 
of J c have been chosen to be equally spaced in logarithm 
scale and they range (from lower to higher maximum force 
—thinner to thicker line) (a) from J c = 1.000 10 6 A/m 2 to 
J c = 5.623 10 7 A/m 2 ; (b) and (c) from J c = 1.000 10 7 A/m 2 to 
J c = 5.623 10 8 A/m 2 . For comparison, the loop corresponding 
to J c = 1.000 10 7 A/m 2 is showed as a dotted line. 



FIG. 6: (a) Vertical stiffness k zz and (b) hysteretic energy 
losses E zz as a function of the aspect ratio of the superconduc- 
tor, for different values of the critical current: J C R/ Ho=0.1 
(solid line), J c R/Ho = l (dot-dashed line). Other parameters 
are written in the text. 



FIG. 7: Force versus distance minor loops started at 
2i=0.005m and of amplitude Az=0. 0025m for the cases 
{a)J c R/H =Q.\ and (b) J C R/H =1. The value of L/R is 
(from lower to higher maximum force — thinner to thicker 
line) L/R = 0.1, 0.5, 1.0, 1.5, and 2.0, respectively. 



FIG. 8: (a) Vertical stiffness k zz and (b) hysteretic energy 
losses E zz as a function of the relative radius of the PM and 
the superconductor a/R, for different values of the critical 
current: J C R/ Hq= 0.1 in dotted line and JcR/Ho — 1 in solid 
line. Other parameters are written in the text. 
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